FIVE TYPES OF BLOW-UP IN A SEMILINEAR 
FOURTH-ORDER REACTION-DIFFUSION EQUATION: 
AN ANALYTIC-NUMERICAL APPROACH 



V.A. GALAKTIONOV 

Abstract. Five types of blow-up patterns that can occur for the 4th-order semilincar 
parabolic equation of reaction-diffusion type 

u t = -A 2 u + \u\ p ~ 1 u in R N x (0,T), p > 1, lim t ^ T - sup^N \u(x, t)\ = +00, 

are discussed. For the semilinear heat equation u t = Au + u p , various blow-up patterns 
were under scrutiny since 1980s, while the case of higher-order diffusion was studied 
much less, regardless a wide range of its application. The types of blow-up include: 

(i) Type I(ss): various patterns of self-similar single point blow-up, including those, for 
which the final time profile \u(-,T')\ N ^ p ~ 1 ^ 4 is a measure; 

(ii) Type I(log): self-similar non-radial blow-up with angular logarithmic TW swirl; 

(iii) Type I(Her): non self-similar blow-up close to stable/centre subspaces of Hermitian 
operators obtained via linearization about constant uniform blow-up pattern; 

(iv) Type Il(sing): non self-similar blow-up on stable/centre manifolds of a singular 
steady state in the supercritical Sobolev range p> ps — jj^i f° r N > 4; and 

(v) Type II(LN): non self-similar blow-up along the manifold of stationary generalized 
Loewner-Nirenberg type explicit solutions in the critical Sobolev case p — ps, when 
|u(«, J 1- JI^Cp - 1 )/ 4 contains a measure as a singular component. 

All proposed types of blow-up are very difficult to justify, so formal analytic and 
numerical methods are key in supporting some theoretical judgements. 



I. From second-order to higher-order blow-up R-D models: a PDE 

ROUTE FROM XXTH TO XXIST CENTURY 

1.1. The RDE— 4 and applications. This paper is devoted to a description of blow-up 
patterns for the fourth-order reaction- diffusion equation (the RDE-4 in short) 

(1.1) u t = -A 2 u+ \u\ p ~ l u in R N x(0,T), where p > 1, 

where A stands for the Laplacian in R . This has the bi-harmonic diffusion —A 2 and is 
a higher-order counterpart of classic second-order PDEs, which we begin our discussion 
with. For applications of such higher-diffusion models, see short surveys and references in 
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[31 HQ]. In general, higher-order semilinear parabolic equations arise in many physical ap- 
plications such as thin film theory, convection-explosion theory, lubrication theory, flame 
and wave propagation (the Kuramoto-Sivashinsky equation and the extended Fisher- 
Kolmogorov equation), phase transition at critical Lifschitz points, bi-stable systems and 
applications to structural mechanics; the effect of fourth-order terms on self-focusing prob- 
lems in nonlinear optics are also well-known in applied and mathematical literature. For 
a systematic treatment of extended KPPF-equations, see Peletier-Troy [76J. 

Note that another related fourth-order one-dimensional semilinear parabolic equation 

(1.2) u t = -u xxxx - [(2 - (u x ) 2 )u x ] x - au + qe su , 

where a, q and s are positive constants obtained from physical parameters, occurs in the 
Semenov-Rayleigh-Benard problem [53], where the equation is derived in studying the 
interaction between natural convection and the explosion of an exothermically-reacting 
fluid confined between two isothermal horizontal plates. This is an evolution equation 
for the temperature fluctuations in the presence of natural convection, wall losses and 
chemistry. It can be considered as a formal combination of the equation derived in [H] 
(see also [5]) for the Rayleigh-Benard problem and of the Semenov-like energy balance 
[79t [T6] showing that natural convection and the explosion mechanism may reinforce each 
other; see more details on physics and mathematics of blow-up in [10] . In a special limit, 
(11.2|) reduces to the generalized Frank-Kamenetskii equation (see [3] for blow-up stuff) 

(1-3) Ut = u xxxx + e , 

which is a natural extension of the classic Frank-Kamenetskii equation; see below. 

Equation fll.ip can be considered as a non-mass-conservative counterpart of the well- 
known limit unstable Cahn-Hilliard equation from phase transition, 

(1.4) u t = -u xxxx - (\u\ p ~ x u) xx in Ex R + , 

which is known to admit various families of blow-up solutions; see [10] for a long list of 
references. Somehow, fll.ll) is related to the famous Kuramoto-Sivashinsky equation from 
flame propagation theory 

(1.5) u t = —u xxxx — u xx + uu x in Kx R + , 
which always admits global solutions, so no blow-up for (11.51) exists. 

1.2. On second-order reaction-diffusion (R— D) equations: a training ground of 
blow-up PDE research in the XXth century. Blow-up phenomena, as examples of 
extremely nonstationary behaviour of nonlinear mechanical and physical systems, become 
more natural in PDE theory since a systematic developing combustion theory in the 1930s. 
This essential combustion influence began with the derivation of the semilinear parabolic 
reaction-diffusion PDE such as the classic Frank-Kamenetskii equation (1938) [15] 

(1.6) u t = Au + e u in R N x R + , 

which occurs in combustion theory of solid fuels and is often also called the solid fuel 
model. First blow-up results in related ODE models are due to Todes in 1933; see the 
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famous monograph [87] for details of the history and applications. The related model with 
a power superlinear source term takes the form (also available among various nonlinear 
combustion models [87J) 

(1.7) ut = Au + \u\ p ~ l u in R N x R + , where p > 1. 

Thus, for such typical models, blow-up means that in the Cauchy proble the classic 
bounded solution u = u(x,t) exists in M. N x (0,T), while 

(1.8) supjjgjjjv \u(x, t)\ — > +oo as t—>T~, 

where T 6 M + = (0, +oo) is then called the blow-up time of the solution u(x, t). 

During last fifty years of very intensive research starting from seminal Fujita results 
in 1966 (on what is now called Fujita exponents) , we have currently got rather complete 
understanding of the types of blow-up for the semilinear (11.61) . (11.71) and other models. 
This is very well explained in a number of monographs; see [2] [751 EH HU (66J [391 121] 177]. 

However, one should remember that even for simple R-D equation such as ( 11.6jl and 
(II. 7ft . there are blow-up scenarios in the multi-dimensional geometries, which still did 
not get a proper rigorous mathematical justification. For instance, there are a number 
surprises even in the radial geometry for (ll.7p . which reads for r = |se| > as 

(1.9) Ut = -phi (r N ~ 1 u r ) r + |m| p_1 w in M + x M + (w r |r=o = 0, symmetry), 
in the supercritical Sobolev range 

(1.10) P>Ps = jt%, where N > 2. 

Several critical exponents, which may essentially change blow-up evolution, appear for 
fll.9p in the range (11.101) . among those let us mention the most amazing ones: 

(1.11) pj L = 1 + N ^J 2 ^- V p L = l + J £ ro (N> 11); p M = 1 + ^ (N > 12); etc. 
In particular, this shows that, in the parameter range 

(1.12) JV>11 and p > p s , 

new principal issues of blow-up evolution for (jl.9p essentially take place. Note that, in 
[38] , some critical blow-up exponents were shown to exist for the quasilinear combustion 
equation with a porous medium diffusion: 

(1.13) u t = Au m + u p , where p > m > 1 (u(x, t)>0). 

This shows certain universality of formation of blow-up singularities for a wider class of 
R-D equations, which now we are going to extend to the RDE-4 (11.11) . 

We do not plan to give any detailed enough review of such a variety of these delicate 
and becoming diverse (rather surprisingly) in the XXIst century mathematical results, 
which quite recently attracted the attention of several remarkable mathematicians from 
various areas of PDE theory. We refer to [83l HHJ [38] for earlier results since 1980s and 
90s, and to more recent papers |T3],[6T] and [67]-[69] as a guide to the research, which was 

1 For simplicity, we avoid using initial-boundary value problems, where boundary conditions can affect 
some manipulations and speculations around; though can be included. 
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essentially intensified last few years. Further results can be traced out by the MathSciNet, 
using most recent papers of the authors mentioned above. 

It is worth mentioning that most of these results have been obtained for nonnegative 
blow-up solutions of (11.6p . (11.71) . and (11.101) . since the positivity property is naturally 
supported by the Maximum Principle (the MP) for such second-order parabolic equations. 
For instance, a full classification of such nonnegative blow-up patterns for (11.71) (all of them 
belong to the family Type I(Her)) in the subcritical range 1 < p < p$ was obtained in 
[83] . For (II. 6p . this happens in dimension N = 1 and 2. In other words, the family 
of blow-up patterns for (11.61) and subcritical (11.71) first formally introduced in [83] is 
evolutionary complete (a notion from [22], where further references can be found). In the 
range p > p$ for (ll.7p and from iV = 3 for (II. 6p . there occur self-similar patterns of Type 
I(ss) and many others being non-self-similar, which makes the global blow-up flow much 
more complicated. 

For p > ps, such a complete classification for (11.71) is far from being complete. E.g., 

(1.14) for p > ps in (II. 7p . nonsymmetric blow-up patterns are practically unknown. 

Moreover, for solutions of changing sign, the results are much more rare and are essen- 
tially incomplete. It is worth mentioning surprising blow-up patterns of changing sign 
constructed in [T4], with the structure to be used later on for (jl.ip . where we comment 
on this Type II(LN) blow-up patterns for (11.11) in greater detail. 



1.3. Back to the RDE— 4: five types of blow-up patterns and layout of the 
paper. We are going to discuss possible types of blow-up behaviour for the RDE-4 (II. ip . 
In what follows, we are using the auxiliary classification from Hamilton [35], where Type 
I blow-up means the solutions satisfying, for some constant C > (depending on u), 

Type I: (T — t)p^\u(x, t)\ <C as t — > T~ , and, otherwise 
(1.15) ^_ 

Type II: lim sup^j.- (T — t)p- x sup x \u[x,t)\ = +oo 

(Type II also called slow blow-up in [3H]). In R-D theory, blow-up with the dimensional 
estimate (11.151) was usually called of self-similar rate, while Type II was referred to as 
fast and non self-similar; see [39] and [78J . 

Thus, we plan to describe the following five types of blow-up with an extra classification 
issues in each of them (this list also shows the overall layout of the paper): 

(i) Type I(ss): various patterns of self-similar single point blow-up mainly in radial 
geometry, including those, for which \u(-, T~)\ N( ^ p ~ lS) ^ is a measure (Section [21 almost 
nothing is known for non-radial similarity blow-up patterns for N > 2, which are expected 
to exist); 

(ii) Type I(log): non radial self-similar blow-up with angular logarithmic travelling wave 
(logTW) swirl, which in the similarity rescaled variables corresponds to periodic orbits 
as aj-limit sets (Section [3]); 
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(iii) Type I(Her): non self-similar blow-up close to stable/centre subspaces of Hermitian 
operators obtained via linearization about constant uniform blow-up and matching with 
a Hamilton- Jacobi region (Section HJ); 

(iv) Type Il(sing): non self-similar blow-up on stable/centre manifolds of singular 
steady state (SSS) in the supercritical Sobolev (and Hardy) range p > ps = for 
N > 4 with matching to a central quasi-stationary region (Section ED; and 

(v) Type II(LN): non self-similar blow-up along the manifold of stationary generalized 
Loewner-Nirenberg type explicit solutions in the critical Sobolev case p = ps, when final 
time profiles \u(-, T~) \ n (p~ 1 )/ a contain measures in the singular component (Section ED- 

We must admit that the analysis of all the blow-up type indicated above is very difficult 
mathematically, so we do not present practically no rigorous results. Recall that, even for 
the second-order equation (11.61) . all these types excluding Type I(Her) still did not have 
not only any complete classification, but some of them were not detected at all. For (11. ip . 
the best known critical exponent is obviously Sobolev's one 

(1.16) ps = |pzf , where N > 5, and also p* = 

while the others, as counterparts of those in (11.111) . need further study and understanding. 
However, many critical exponents for (II. ip cannot be explicitly calculated. Overall, we 
aim that our approaches to blow-up patterns can be extended to 2mth-order parabolic 
equations such as 

(1.17) u t = -(-A) m u+\u\ p - 1 u, 

though the case m = 2 (the first even m's) already contains some surprises. 

Nevertheless, it seems that, at some stage of struggling for developing new concepts, 
it is inevitable to attempt to perform a formal classification under the clear danger of a 
lack of any rigorous justification^. In this rather paradoxical connection, it is also worth 
mentioning that the most well-known nowadays and the fundamental open problem of 
fluid mechanics^ and PDE theory on global existence or nonexistence (blow-up) of bounded 
smooth L 2 -solutions of the Navier-Stokes equations (the NSEs) 

(1.18) u t + (u- V)u = -Vp + Au, divu = in l 3 xl + , u| t=0 = u G L 2 n 

from one side belongs to a "blow-up configurational" type: to predict possible swirling 
"twistor-tornado" type of blow-up patterns. Moreover, it seems that the NSEs (I1.18P was 
the first model, for which J. Leray in 1934 [571 P- 245] formulated the so-called Leray's 
scenario of self-similar blow-up as £ ^ T~ and a similarity continuation beyond for 
t > T. Nonexistence of such similarity blow-up for the NSEs (I1.18P was proved in 1996 in 
Necas-Ruzicka-Sverak [71]. However, for the semilinear heat equations (11.71) and fl!.13p . 
the validity of Leray's scenario of blow-up was rigorously established; see [38, 68J and 
references therein. 

2 Actually following Kolmogorov's legacy from the 1980s sounding not completely literally as: "The 
main goal of a mathematician is not proving a theorem, but an effective investigation of the problem..." . 
3 The Millennium Prize Problem for the Clay Institute; see Fefferman [T^] . 
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In general, we observe certain similarities between these two blow-up problems; see 
[2"T] . where Type I (log) patterns were introduced for (ll.ISp and [2S] for more details and 
references on other related exact blow-up solutions. Overall, we claim that equations fll .6f> . 
fll.7p . (11. II) . and (11.181) admit some similar principles of constructing various families of 
blow-up patterns, though, of course, for the last two ones, the construction gets essentially 
harder and many steps are made formally, without proper justification. Especially for the 
NSEs (I1.18p . which compose a nonlocal solenoidal parabolic equation: 

(1.19) u ( + P(u-V)u = Au, 

where the integral operator F = I — VA _1 (V-) is Leray-Hopf's projector onto the 
solenoidal vector field. More precisely, the RDE-2 ( II .71) obeying the MP is indeed too 
simple to mimic any NSEs blow-up patterns, while (11. ip . which similar to (I1.19P traces 
no MPs, can be about right (possibly, still illusionary) . Then (II. ip stands for an auxiliary 
"training ground" to approach understanding of mysterious and hypothetical blow-up for 

(HHD. 

In Appendix A, we present other families of PDEs, which expose a similar open problem 
on existence/no nexistence of L°°-blow-up of solutions from bounded smooth initial data. 
Overall, it is worth saying that the problem of description of blow-up patterns and their 
evolution completeness takes and shapes certain universality features in general PDE 
theory of the twenty first century. 



2. Type I(ss): SELF-SIMILAR BLOW-UP 

This is the simplest and most natural type of blow-up for scaling invariant equations 
such as (II. ip . where the behaviour as t — > T~ is given by a self-similar solution: 

(2-1) u s (x,t) = (T-ty^f(y), y= w ^r 7 i, 

where a non-constant function / ^ is a proper solution of the elliptic problem: 

(2.2) A (/)^-A 2 /-h-V/-^/+|/r 1 / = in R N , /(oo) = 0. 

We recall that, for (11.71) . such nontrivial self-similar Type I blow-up is nonexistent in the 
subcritical range p < But this is not the case for the RDE-4 ( 11. ip . Note that ( 12.21) is 
a very difficult elliptic equation with the non-coercive and non-monotone operators, which 
are not variational in any weighted L 2 -spaces. There are no still any sufficiently general 
results of solvability of ( 12. 21) in higher dimensions, so our research is a first attempt. 

In what follows, for any dimension iV > 1, by fo{y) we will denote the first monotone 
radially symmetric blow-up profile, which, being on the lower iV-branch (so fo is not 
unique, see explanations below) is expected to be generic (i.e., structurally stable in the 
rescaled sense). We also deal with the second symmetric profile fi{y), which seems to 
be unstable, or, at least, less stable than f . There are also other similarity solutions 
concentrated about the singular SSS U(y) (see Section [5]), but those, being adjacent to 
the unstable equilibrium U are expected to be unstable also. 
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It is worth mentioning that self-similar blow-up for ( 11. ip is incomplete, i.e., blow-up 
solutions, in general, admit global extensions for t > T. Such principal questions are 
studied in [29] and will not be treated here. 

2.1. One dimension: first examples of nonuniqueness. Thus, for N = 1, ( 12.21) 

becomes the ODE 

(2.3) A(f) = -fW-\yf--± i f + \f\'- l f = in R, /(±oo) = 0, 

which was studied in [5] by a number of analytic-branching and numerical methods. It 
was shown that H 2 . 3 1) admits at least two different blow-up profiles with an algebraic 
decay at infinity. See [4*0~1 § 3] for further centre manifold-type arguments supporting this 
multiplicity result in a similar 4th-order blow-up problem. Without going into detail of 
such a study, we present a few illustrations only and will address the essential dependence 
of similarity profiles f(y) on p. In Figured! we present those pairs of solutions of (12.31) for 
p — | and p = 2. All the profiles are symmetric (even), so satisfy the symmetry condition 

(2.4) /'(0) = f" (0) = 0. 

No non-symmetric blow-up was detected in numerical experiments (though there is no 
proof that such ones are nonexistent: recall that "moving plane" and Aleksandrov's Re- 
flection Principle methods do not apply to ( 11.11) without the MP). Figure [2] shows similar 
two blow-up profiles for p = 5. 

2.2. On existence of similarity profiles for N = 1: classification of blow-up and 
oscillatory bundles. We now provide extra details concerning existence of at least a 
single blow-up profile f(y) satisfying (12.31) . (12.41) . We perform shooting from y = +oo by 
using the 2D bundle (I2.18P to y = 0, where the symmetry condition (12. 4p are posed (or to 
y = — oo, where the same bundle ( I2.18P with y \— > —y takes place). By / = f(y; C\, C2), 
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u t =-u xxxx + ' U ' P tw0 blow-up similarity profiles for p=5 




y 



FIGURE 2. Two self-similar blow-up solutions of (12. 3D for p = 5. 

we denote the corresponding solution defined on some maximal interval 

(2.5) y G (y , +oo), where y = Vo{Ci, C 2 ) > -oo. 

If yo(Ci) C 2 ) = —oo, then the corresponding solution f(y; C\, C 2 ) is global and can repre- 
sents a proper blow-up profile (but not often, see below). Otherwise: 

(2.6) y (C 1 ,C 2 ) > -oo f(y; C u C 2 ) -»• oo as y -»• j/J. 
Note that "oscillatory blow-up" for the ODE close to y = t/g : 

(2.7) / W = |/r 1 /(l + o(l)), 

where limsup f(y) = +oo and liminf f(y) = —oo as y —>■ y$, is nonexistent. The proof is 
easy and follows by multiplying (12. 7p by /' and integrating between two extremum points 
(yi, 2/2), where the former one y% is chosen to be sufficiently close to the blow-up value 
whence the contradiction: 

o<Un 2 (yi)--^- 1 \f\ p+1 (yi)<o. 

We first study this set of blow-up solutions. These results are well understood for such 
fourth-order ODEs; see [12] , so we omit some details. 

Proposition 2.1. The set of blow-up solutions (12. 7p is four- dimensional. 

Proof. The first parameter is yo G R. Other are obtained from the principal part of the 
equation (12. 7\i describing blow-up via (12.31) as y — >• y$. We apply a standard perturbation 
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Funsction 4>(m)=m(m-1)(m-1)(m-3) 
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FIGURE 3. The graph of function <3?(m) in (|2.8|) . and towards solutions of (|2.10p . 

argument to (I2.7I) . Omitting the o(l)-term and assuming that / > 0, we find its explicit 
solution 

(2.8) f (y) = A (y - y )-&, A*' 1 = $(-^), $(m) = m(m - l)(m - 2)(m - 3). 

For convenience, the graph of $(m) is shown in Figure [31 Note that it is symmetric 
relative to m = |, at which $(m) has a local maximum: 

(2.9) $(§) = A. 
By linearization, / = Jo + ^ we get Euler's ODE: 

(2.10) (y - 2/o) 4 r (4) = pAl 1 = paK-jjij). 

It follows that the general solution is composed from the polynomial ones with the fol- 
lowing characteristic equation: 

(2.11) Y(y) = (y-y )™ => $(m) = ^(-^). 

Since the multiplier p > 1 in the last term in (12.1 ip and m = — ^-j- is a solution if this 
"px" is omitted, this algebraic equation for m admits a unique positive solution m + > 3, 
a negative one m_ < — ^y, which is not acceptable by (I2.8p . ad two complex roots my 

with Re my = | > 0. Therefore, the general solution of (12. 7p about the blow-up one 
( 12.81) . for any fixed ?/o, has a 3D stable manifold. □ 

Thus, according to Proposition 12.11 the blow-up behaviour with a fixed sign ( \2.6h (i.e., 
non-oscillatory) is generic for the ODE (12.31) . However, this 4D blow-up bundle together 
with the 2D bundle of good solutions (I2.18P as y — > ±oo are not enough to justify the 
shooting procedure. Indeed, by a straightforward dimensional estimate, an extra bundle 
at infinity is missing. 
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FIGURE 4. A periodic solution of = I/]* -1 / separatrix: p = 5. 



To introduce this new oscillatory bundle, we begin with the simpler ODE (I2.7p . without 
the o(l)-term, and present in Figure H] the results of shooting of a "separatrix" that lies 
between orbits, which blow-up to ±oo. Obviously, this separatrix is a periodic solution of 
this equation with a potential operator. Such variational problems are known to admit 
periodic solutions of arbitrary period. 

Thus, Figure H] fixed a bounded oscillatory (periodic) solution as y — > +oo. When we 
return to the original equation (I2.3I) . which is not variational, we still are able to detect 
a more complicated oscillatory structures at y — oo. Namely, these are generated by the 
principal terms in 

(2.12) / (4) = _i fy+ | / | P -i / + ... as y^cc. 

Similar to Figure HI in Figure [5], we present the result of shooting (from y = — oo, which 
is the same by symmetry) of such oscillatory solutions of (12.31) for p — 5. It is easy to see 
that such oscillatory solutions have increasing amplitude of their oscillations as y — > oo, 
which, as above, is proved by multiplying ( 12.121) by /' and integrating over any interval 
2/2) between two extrema. Figure [6] shows shooting of similar oscillatory structures at 
infinity for p = 7 (a) and p = 2 (b). It is not very difficult to prove that the set of such 
oscillatory orbits at infinity is ID and this well corresponds to the periodic one in Figure 
H] depending on the single parameter being its arbitrary period. 

By Cj (Ci) in Figure O we denote the values of the second parameters C 2 such that, 
for a fixed C\ G 1R, the solutions f(y; Ci, Cif ) blow up to ±00 respectively. These values 
are necessary for shooting the symmetry conditions (12. 4p . 

10 




Shooting oscillatory solution: p=5, N=1 , C =1 , C = 2.8956362303297 
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FIGURE 5. Shooting an oscillatory solution at infinity of (|2,3|) : p = 5. 



Schooting oscillatory solution: p=7, N=1, (^=1, C 2 =2. 941 16376644039 



Schooting oscillatory solution: p=2, N=1, ^=1, C 2 =0. 1131218908 
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(a) p = 7 



(b) p = 2 



FIGURE 6. Shooting an oscillatory solution at infinity of (|2.3|) : p = 7 (a) and 
p = 2 (b). 



(2.13) 



find Ci 2 such that: 



Thus, overall, using two parameters Ci j2 in the bundle (12.181) for y ^> 1 leads to a 
well-posed problem of a 2D-2D shooting: 

yo(C 1 ,C 2 ) = -oo, and 
no oscillatory behaviour as y — >• — oo. 

Concerning the actual proof of existence via shooting of at least a single blow-up pat- 
terns fo(y), by construction and oscillatory property of the equation (12.31) . we first claim 
that in view of continuity relative to the parameters, 

(2.14) for any C x > 0, there exists C 2 *(Ci) G (C 2 (d), C 2 + (Ci)) such that f"'(0) = 0. 

n 



We next change C\ to prove that at this C|(Ci) the derivative /'(0) also changes sign. 
Indeed, one can see that 

(2.15) /'(0;Ci,C*(Ci)) > for d < 1 and /'(0; Ci, C|(Ci)) < for Ci > 1. 

Actually, this means for such essentially different values of C\, the solution f(y; C\, C%(Ci) 
has first oscillatory "humps" for y > and y < respectively. By continuity in C\, (12.151) 
implies existence of a C* such that 

(2.16) f(0;C*,C 2 *(C*)) = 0, 

which together with (I2.14p induced the desired solution. Overall, the above geometric 
shooting well corresponds to that applied in the standard framework of classic ODE 
theory, so we do not treat this in greater detail. However, we must admit that proving 
analogously existence of the second solution fo(y) (detected earlier by not fully justified 
arguments of homotopy and branching theory and confirmed numerically) is an open 
problem. A more difficult open problem is to show why the problem (12.131) does not 
admit non-symmetric (non-even) solutions f(y) (or does it?). 

2.3. Dimensions N > 2: on 2D shooting and analogous nonuniqueness. In higher 
dimensions, it is easier to describe Type I(ss) blow-up in radial geometry, where (12.21) also 
becomes an ODE of the form (now y stands for \y\ > 0) 

(2.17) -/« - ^ f >» - (n^w^i r + (*-y-3) f _± yf __L_ f+ lfr i f = o 

in R, with the same two symmetry condition (12.41) . To explain the nature of difficulties 
in proving existence of solutions of ( 12. XTD . let us describe the admissible behaviour for 
j/ > 1. There exists a 2D bundle of such asymptotics (see details in [TT], § 3.3]): as 
y -> +oo, 

(2.18) f(y) = [C x y-^ + ...] + [C 2 y~^ N ~^ ] e - a ^ 4/3 + ...], where a = 3 • 2~§ 

and C\ and C 2 are arbitrary parameters. This somehow reminds a typical centre manifold 
structure of the origin {/ = 0} at y = 00: the first term in ( 12. 18ft is a node bundle with 
algebraic decay, while the second one corresponds to "non-analytic" exponential bundle 
around any of algebraic curves. Thus, a dimensionally well-posed shooting is: 



(2.19) Shooting: using 2 parameters Ci t2 in (1 2 . 1 8 [) to satisfy 2 conditions (12. 4p . 



In case of analytic dependence of solutions of ( 12.171) on parameters in the bundle 
(12.181) (this is rather plausible via standard trends of ODE theory, but difficult to prove), 
the problem cannot have more than a countable set of solutions. Actually, our numerics 
confirm that in wide parameter ranges of p > 1 and N > 1, there exist not more than 
two solutions (up to other more unstable ones about the SSS; see Section [5]): 

(2.20) fo(y) with {Cx (p,N),C 2 o(p,N)}, and f x (y) with {^(p, N),C 21 (p, N)}. 



The rest of this section is devoted to justify this. 
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u =-A 2 u+|u| p 'u: two blow-up similarity profiles for p=2 and N=3 u f =-A 2 u+ |u| p 1 u: two similarity profiles for p=2 and N=10 




y y 

(a) p = 2, N = 3 (b) p = 2, N = 10 



FIGURE 7. Two self-similar blow-up solutions of (12.17P for p = 2: N = 3 (a) 
and JV = 10 (b). 



The eventual similarity blow-up patterns can be characterized by their final time pro- 
files: passing to the limit t — > T~ in ( 12. ip and using the expansion ( 12. 18ft yields 

(2.21) if Ci(p,N) ^ 0, then u(z, i) -> Ci|x|"^ as £ _, t~ 

uniformly on any compact subset of R \ {0}. If Ci = in (12.181) . i.e., f(y) has an 
exponential decay at infinity, then the limit is different: in the sense of distributions, 



(2.22) d(p,iV) = 



* -> C 3 «y(x), t -> T ; C 3 = /|/| * < oo. 



It is very difficult to prove that (12.221) actually takes place at some p = ps{N) > 1 (even for 
N — 1), and we will justify this numerically for some not that large dimensions N < 11. 

We now start describe various similarity blow-up profiles for N > 2. As a first and 
analogous to N — 1 example, in Figure [TJ we construct numerically first two profiles, 
fo(y) and fi{y), for the three-dimensional case N = 3 and iV = 10 for p = 2, which look 
rather similar to those in Figure [1] for N — 1 . 

2.4. > 2: p-branches of the profile fo(y) and fi(y). Such p-branches of solutions 
are a convenient way to describe families of profiles fo{y) depending on the exponent p; 
cf. [3D] |4"T] . In Figure El we present such a branch of /o for = 4, where (a) shows 
the actual smooth deformation of fo(y) with changing p, while (b) is the corresponding 
p-branch. In Figure O the same is done for N = 8. Note that both Figures (b) show that 
\\f\\oo = /(0) approaches 1 for large p, which is a general phenomenon for such ODEs 
described in [30, § 5]. Similarly, Figure fTOl shows p-branches of the second blow-up profile 
fi{y) for p = 2 in the cases N = 1 (a) and iV = 12 (b) (the critical dimension, where 
Ps — 2). 
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f(y) 



p-branch for N=4: deformation of f Q (y) for p g [2,25] 



p-branch of f (y) for N=4 and p e [2,25] 




10 15 
P 



(a) /o-deformation (b) p-branch 

FIGURE 8. The p-branch of blow-up self-similar profiles fo(y) for N = 4. 



p-branch for N=8: deformation of f (y) for p e [2,45] 




p-branch of f for N=8: p e [2,45] 



5 10 15 20 25 30 35 40 45 

P 



(a) /o(y)-defoi'mation 



(b) p-branch 



FIGURE 9. The p-branch of blow-up self-similar profiles fo(y) for N = 8. 



It is well understood that for equations such as (I2.17p . the solutions f(y) blow-up as 
p — > 1 + with a super-exponential rate ~ (p — 1)~ 1 /(p _1 ); see [501 |4"T] . As an example, in 
Figure dU we present such a blowing up behaviour of the p-branch of fi{y) for N = 6. 

2.5. N > 2: iV-branches of blow-up profiles. Firstly, in some iV-intervals, there is a 
continuous dependence of fo{y) on the dimension, as Figure [T21 clearly shows for p = 2 
and Figure [12] for p = 5 (those values of p will be constantly used later on for the sake of 
comparison) . 

However, we found that there are other solutions of the monotone type fo, which are 
shown in Figure [TH for p = 5 (a) and p = 2 (b), where in the latter one the profile from 
the lower iV-branch is not shown as being too relatively small. 
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p-branch of f, (y) for N=1 : p e [2,20] 



p-branch of f (y) for N=12: p E [2,20] 



10 
P 



(a) p-branch for N = 1 




(b) p-branch for N = 12 



FIGURE 10. The p-branches of blow-up self-similar profiles fi(y) for p = 2: 
N = 1 (a) and N = 12 (b). 



Blow-up of p-branch of f (y) as p tends to 1 + : N=6 



Blow-up of p-branch of f (y) as p tends to 1 + for N=6 



f(y) 




(a) /i-deformation 




1.2 1.3 1.4 1.5 1.6 1.7 1 



(b) p-branch 



FIGURE 11. Blow-up as p — > 1 of the p-branches of profiles fi(y) for N = 6: 
p-deformation of fi(y) (a) and blow-up of p-branch in the log-scale (b). 



Thus, secondly, this nonuniqueness demands another approach to branching, namely, 
the iV-branching that we perform next. In Figure [151 we show the lower iV-branch of 
solutions fo(y) for p = 5, where (a) describes the deformation of f and (b) gives the 
actual p-branch. Blow-up of the upper iV-branch as 

(2.23) N^N+, where N p : $±f = p, 

so that A^5 = 6 for p = 5, is shown in Figure [T6| with the same meaning of (a) and (b). A 
general view of the whole iV-branch of fo for p = 5 is schematically explained in Figure 
rm where by dotted line we draw a possible expected but still hypothetical connection of 
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u ( =-A 2 u+ |u| p 1 u: first similarity profile f Q (y) for p=2 and N=1,2,...,16 




2 4 6 8 10 12 



y 

FIGURE 12. Similarity blow-up solutions fo(y) of (12.170 for p = 2 in dimensions 
N= 1, 2, 16. 




y 

FIGURE 13. Similarity solutions f of OT?]) for p = 5 and N = 1, 2, ... , 15. 

the lower (stable) and the upper (more unstable, plausibly) / -branches, which we were 
not able to reconstruct numerically. Numerical continuation in the parameter iV is quite 
a challenging problem in some iV-ranges. 
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Two blow-up profiles for p=5, N=7 



Monotone blow-up profile f (y) from upper N-branch: p=2, N=13 




8 10 



(a) p = 5 




(b) p = 2 



FIGURE 14. Two monotone blow-up profiles of /o-type: for p = 5 (a) and for 
P = 2(b). 



N-branch for p=5: smooth deformation of f (y) for n e[1 ,14.4] 




N-branch of f (y) for p=5 for Ne [1 ,14.4] 



(a) /o-deformation 



(b) lower iV-branch 



FIGURE 15. Lower iV-branch of f for p = 5. 



Thus we expect that there exists a saddle-node bifurcation at some 

(2.24) p = 5: N sn e (14.979,15). 

In Figure fT8l for p = 2, we show blow-up of the upper A^-branch as iV — >• A^2 = 12 + . We 
then expect that the lower and upper branches have a turning (saddle-node) bifurcation 
point at some 

p = 2: N sn G (20.3, 23]. 

We hope that such an interesting saddle-node branching phenomenon will attract true 
experts in numerical methods, bearing in mind that numerical experiments might be for 
a long time the only tool of the study of such blow-up phenomena. 
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0.5 1 1.5 2 6 6.2 6.4 6.6 6.8 7 

y N 

(a) /o-deformation (b) upper ./V-branch 



FIGURE 16. Upper iV-branch of /o for p = 5 with blow-up as N — > 6 + . 




FIGURE 17. The whole iV-branch of / for p = 5. 



Finally, in Figure [T9], we present the numerical results confirming iV-branching for 
p = 2 of the second blow-up profile fi, where, as usual, (a) describes smooth deformation 
of fi(y), while (b) shows the iV-branch. It seems that iV-branches of fi are global and do 
not suffer from a saddle-node bifurcations. 

2.6. On sign changes of fo(y) and fi(y). We now study some particularly important 
properties of blow-up similarity profiles. We begin with the easier property of sign changes. 
We have seen already several strictly positive profiles fo(y) for some p's, which is rather 
surprising since the equations do not obey the Maximum Principle. However, we will 
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N-branch for p=2 with N e [1 ,20.3] 



Blow-up of upper N-branch for p=2 near N =12 + 




(a) lower TV-branch 
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N 



(b) upper TV-branch 



FIGURE 18. Lower and upper TV-branches of /o for p = 2: the lower branch for 
TV S [1, 20.3] (a) and blow-up of the upper one as TV — > 12 + (b). 



N-branch for p=2: deformation of f (y) for N e [1 ,20] 




N-branch of f (y) for p=2: Ne [1,20] 




(a) /o-deformation (b) TV-branch 

FIGURE 19. The TV-branch of fi for p = 2 and TV G [1, 20]. 

show that, for smaller p, the similarity profiles can gain extra zeros as sign changes. Such 
p's , when a zero is gained, we denote by p (N). 

Consider one dimension TV = 1. Firstly, the attentive Reader can see that in Figure 
COX), already for p — | the profile fo(y) changes sign, while for p — 2, it is positive. Hence, 
Po{l) G (|,2). Secondly, more thorough numerics are presented in Figure EH where (a) 
shows Jo's in a vicinity of 

(2.25) po(l) = 1.7358... , 

while (b) shows a sharp shooting of the critical value f 1 2 . 2 5 j) . Figure EH shows shooting 

rf ) (l) = 1.23..., 
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at which the second profile f\{y) gets a new zero. By the boldface line we denote a new 
"p-undetected" solution with extra zeros gained at another p i> showing that such roots 
are not unique. Since f\(y) is expected to be less stable, we will concentrate on the roots 
Po for the generic blow-up profile fo(y)- 
Thus, similarly, Figure 1227 a) yields 

(2.26) p (3) = 1.446... and p (9) = 1.204... . 

In Figure 1237 a). a similar phenomenon is checked for iV = 10, with po(10) = 1.188... . In 
(b), we see no sign changes of fo(y) for N = 12, but this happens for smaller p « 1 + , when 
/(0) gets 10 5 -10 6 , while their negative counterparts take values ~ — 10 5 , and numerics 
become rather unreliable. The overall numerical results for shooting po(N) are shown in 
Table 1. 

Table 1. Roots po{N) for profiles fo(y) to get a new zero 



N 


Po(N) 


1 


1.7358... 


2 


1.53... 


3 


1.446... 


4 


1.377... 


5 


1.320... 


6 


1.28... 


7 


1.25... 


8 


1.226... 


9 


1.204... 


10 


1.188... 


11 


1.16... 



A proper asymptotic theory for p m 1 + involving expansions such as (I2.18P and (p — 
l)"^r-scaling of the ODE lETTj) (cf. [H § 5]) would be fruitful. Note that we have 
observed some numerical evidence for existence of the second root poi f° r N = 8 and 9 
(see the dotted line in Figure 122!) . which is surprising in view of non-oscillating of the 
exponential term in (I2.18p . but numerics were too difficult and rather poor to identify the 
new root if any. 

2.7. On final time measure-like Type I(ss) blow-up. This is a much more difficult 
problem, which we resolve numerically for N = 1 only. For iV > 2, we got no sufficiently 
reliable results (rather plausibly, such a self-similar phenomenon may be unavailable in 
some higher dimensions). 

We will refine Figure I2"U1 We claim that the equation (I2.22p has the following root: 

(2.27) p s (l) = 1.40... < p (l) = 1-7358... , 

at which the coefficient Ci(p, N) vanishes, so that fo(y) has exponential decay at infinity. 
To see this, we show in Figure [241 with the scale of 10 -20 how the coefficient C\{jp, 1) 
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(a) /o profiles, N = 1 (b) shooting po{l), N = 1 

FIGURE 20. Shooting the root p Q (l) for N = 1. 



N=1, profile f (y): p =1.155 p' 1| =1.23... 



_ 3 N=1, profile f^y): p mi =1.155 p^ >=1 .23... 





(a) h profiles, N = 1 (b) shooting p^l), N = 1 

FIGURE 21. Shooting the root ^(l) for N = 1. 



changes sign around (12. 2 7ft : 

(2.28) Ci(1.39, 1) > 0, while d(1.40, 1) < 0. 

Non vanishing of these two profiles in smaller scaled up to 1CT 40 was checked in the 
logarithmic scale (we do not present here a number of such numerics). 

2.8. On non-radial self-similar blow-up patterns in dimensions N > 2. This 
question was not studied in the literature at all and indeed is very difficult. We make 
a slight observation only: the performed below linearization (14. 2 j) about the constant 
equilibrium in the elliptic equation (12. 21) leads to the perturbed linear elliptic equation 



(2.29) 



(B* + I)Y + D(Y) = 
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_ 4 N=3: shooting p (3)=1. 446... of sign change of f (y) 




N=9: shooting p (9)=1 .204..., f(0)- 12000 




(a) shooting po(3) 



(b) N = 9 



FIGURE 22. Shooting the root po(l)- for N = 3 (a) and for N = 9 (b). 



_5 N=10: shooting p (10)=1.188..., f(o) -62000 




N=12: f (y)>0 for p<2, "nonexistence" of p n (12) 




(a) shooting po for N = 10 



(b) N = 12 



FIGURE 23. Shooting the root po(l): for iV = 10 (a) and no sign changes of 
f (y) for N = 12 for p > 1.3 (b). 



(on spectral properties of B*, see Lemma l47Tj) . Then, B* + J has a large unstable subspace 
(2.30) E u (0) = Span-O/3 : + 1 > 0}, 

so that the corresponding eigenfunctions may characterize possible shapes of various sim- 
ilarity solutions (actually, this is true for N = 1 [3\). Roughly speaking, we claim that: 



(2.31) 



the dimension M(p, N) = dimi? u (0) — (N + 1) can characterize 
the total number of blow-up similarity patterns as solutions of (12.21) . 



Note that we subtract (N + l)-dimensions corresponding to natural instabilities relative 
to shifting the blow-up point G M> N (N dimensions) and blow-up time T (1 dimension). 
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N=1 : shooting p g =1 .40... such that C^p.N^O 
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y 

FIGURE 24. Shooting the root (j2.27p of the algebraic equation (|2.22j) for N = 1. 

These unstable modes are not available if the blow-up point (0, T) is fixed. The dimension 
M = M(p, N) can characterize the total number of solutions f(y) of the elliptic problem 
(12. 2p including many non-radial ones. In other words, we expect that those unstable 
M modes initiate heteroclinic connections through the corresponding unstable manifold 
Vy u (0) to the set of steady solutions {fk(y), k — 1,2, ...,M}. In Section [U we show that 
stable modes from E s (0) with \p + 1 > and the centre ones E c (0) with Xp + 1 = will 
lead to other "linearized" blow-up patterns, so that {fk} are "nonlinear eigenf unctions" . 

Proving any part of the claim ( 12. 31ft is a difficult open problem for any N > 2. Note 
also that, for the second-order quasilinear counterpart (I1.13P (m > 1 is essential!), non 
radially symmetric self-similar blow-up patterns have been known for more than thirty 
years; see [51] and a survey [SS] for extra details. 

3. Type I(log): SELF-SIMILAR PATTERNS WITH ANGULAR LOGTW SWIRL 

This is a simple idea for producing non-radial blow-up patterns, but its consistency is 
quite questionable. 

3.1. Nonstationary rescaling. Dealing with non-self-similar blow-up, instead of (12. ip . 
we use the full similarity scaling: 

(3.1) u{x,t) = {T-t)-^v(y,T), y = w ^ T7I , r = - ln(T - t) -> +oo, t -> T~ . 
Then v(y,r) solves the following parabolic equation: 

(3.2) v T = A{v) = -A 2 v - \y -Vv - +\v\ p ~ 1 v in R^x^oo), r = -lnT, 
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where A is the stationary elliptic operator in (I2.2p . so that similarity profiles (if any) are 
just stationary solutions of (13.21) . 



3.2. Blow-up angular swirling mechanism. We begin with N = 2, where y = (y±, y?), 
and, in the corresponding polar coordinates {p, tp}, with p 2 = y\ + yf, 

(3.3) A = A P + ^D 2 , where A p = D 2 + \ D p and y • V = pD p . 

We next consider a TW in the angular direction by fixing the angular dependence 

(3.4) ip = ar + p = -a\n(T-t) + p, pe(0,2n), 

where a G 1R is a constant (a nonlinear eigenvalue). In the original independent variables 
{x, t}, (13 .4p represents a blowing up logarithmic TW in the angular direction with unknown 
wave speeds a. In other words, H3.4[) assumes that blowing up as t —>■ T~ is accompanied 
by a focusing TW-angular behaviour also in a logarithmic blow-up manner. 

Thus, assuming the logTW angular dependence (13.41) of the solution v = v(y,p,r), 
ip = ot + p, yields the equation 

(3.5) v T = A(v) - o-Vft = -A 2 v - \ y ■ Vv - ^ v + \v\ p ~ 1 v - av^ in R N x (r , oo), 

where To = — In T. In particular, this non-radial self-similar blow-up may be generated 
by bounded steady profiles satisfying 

(3.6) A(/) - a/, = - A V - 1 y . V/ - / + | / ^ V - <t/„ = in M 2 . 

For a ytz 0, which, as we have mentioned, plays the role of a nonlinear eigenvalue, the blow- 
up behaviour with swirl corresponds to periodic orbits as cj-limit sets; see a discussion 
in [27] to the NSEs (11.181) . As a first approach to solvability of (13.61) . one can assume 
branching of a solution f(p,p) from the radial one fo at a = 0. Then setting / = 
/o + o~ifj* + ... yields that "0*(p, p) must be a nontrivial non-radial eigenfunction for A = 0: 

(3.7) A'(f Q )r = 0. 

On the other hand, branches of solutions / of (13.61) may occur at a saddle-node bifurcation 
a = o"* 7^ 0, where a* belongs to spectrum of the linear pencil A'(f) — aD p . Both 
eigenvalue problem are very difficult, and we do not exclude the possibility that, overall, 
the problem (13. 6p for o ^ may admit such solutions only that are singular at the origin 
y = 0. Anyway, even in this unfortunate case, we believe that introducing such rather 
unknown types of non-radial blow-up with swirl deserves mentioning among other more 
practical patterns. Let us also mention that, in M , one can distribute the variables as 

y=(yuy2,y')eR N , where y E R N ~ 2 , 

and arrange a ovlogTW in variables (2/1,2/2) on ly to get periodic blow-up behaviour. 
Choosing other disjoint pairs (yk,yk+i) and constructing the corresponding periodic swirl 
in these variables, in particular, it is formally possible to produce a quasi-periodic blow-up 
swirl with arbitrary number a\,...,a n , n < [y] , of fundamental frequencies. Of course, 
this leads to complicated nonlinear eigenvalue problems, which are open even for n = 1, 
i.e., for the periodic motion introduced above first. 
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3.3. Remark: on the origin of logTWs and invariant solutions. The scaling group- 
invariant nature of such logTWs seems was first obtained by Ovsiannikov in 1959 |73j . 
who performed a full group classification of the nonlinear heat equation 

u t = (k(u)u x ) x , 

for arbitrary functions k(u). In particular, such invariant solutions appear for the porous 
medium and fast diffusion equations for k(u) = u n , n 7^ 0: 

u t = (u n u x ) x =S> 3 u(x,t)=t^f(x + a\iit), where - ± / + of = (f»f)'. 

Blow-up angular dependence as t — > T~ such as in (I3.4p was studied later on in [1] , where 
the corresponding similarity solutions for the reaction-diffusion equation with source 
(11.131) in M 2 x (0, T) were indicated by reducing the PDE to a quasilinear elliptic problem 
(it seems, there is no still a rigorous proof of existence of such patterns). For parabolic 
models such as (11.13)) . that are order-preserving via the MP and do not have a natural 
"vorticity" mechanism, such "spiral waves" as t —>■ T~ must be generated by large enough 
initial data specially "rotationally" distributed in M 2 . For the biharmonic operator as in 
(II. ip with no MP, such a swirl blow-up dependence may be more relevant; see below. 



4. Type I(Her): NON SELF-SIMILAR "LINEARIZED" PATTERNS WITH A LOCAL 
GENERALIZED HERMITE POLYNOMIAL STRUCTURE 

For the classic R-D equation (|1.7p . a countable set of non self-similar blow-up patterns 
of a similar structure was first formally introduced in [M] , though the history of such non- 
self-similar blow-up asymptotics goes back to Hocking-Stuartson-Stuart in 1972, [50J, who 
invented an interesting novel formal technique of analytic expansions (in fact, an analogy 
of a centre manifold analysis) to confirm that blow-up occurs on subsets governed by the 
"hot spot" variables, as t — > T~: 

(4.1) u(x,t) = (T-ty^[f*(0 + o(l)], where g = -^=^==^ 

and /* is a unique solutions of a Hamilton- Jacobi equation of the form (14.151) (with 
1 1— > ~). A justified construction of such patterns and other applications were performed 
a few years later mainly in dozens of papers by Herrero and Velazquez; see pF?, 83J as a 
guide together with other papers traced by the MathSciNet. It is curious that earlier, in 
1987, a sharp upper bound of such a non-similarity blow-up evolution (14.11) ("first half of 
blow-up") was proved in [SB] by a modification of Friedman-McLeod gradient estimate 
[T7] , though the "second half of blow-up" took extra ten years co complete along similar 
lines, HB1 § 7]. 

For the RDE-4, there is no hope to get an easy and fast rigorous justification of such 
non-self-similar blow-up scenario, though the main idea remains the same. We follow [19] 
and also [20], where such a construction applied to non-singular absorption phenomena 
(regular flows with no blow-up), so a full mathematical justification is available therein. 
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4.1. Linearization and spectral properties. The construction of such blow-up pat- 
terns is as follows. Performing the standard linearization about the constant equilibrium 
in the equation ( 13. 2ft yields the following perturbed equation: 

(4.2) v = /* + Y, /* = (p - => F T = (B* + J)Y + D(y), 

where D(Y) = coY 2 + Co = |(p — l)^ 1 , is a quadratic perturbation as Y — > and 

(4.3) B* = -A 2 -±yV in p%) = e -*\y\ A/ \ oe(0,3-2-f), 
is the adjoint Hermite operator with some good spectral properties j9]: 

Lemma 4.1. B* : i?r,(R ) — > L 2 *(IR Ar ) a bounded linear operator with the spectrum 

(4.4) or(B*) = {A^ = -Jfl , |/3| = 0,1,2,...} ( = <r(B), B = -A 2 + \ y ■ V + f /). 
Eigenfunctions ij)%(y) are \(3\th-order generalized Hermite polynomials: 

(4.5) ^(y) = My" + E?fi /4] J(A) V] , |/?| = 0, 1, 2, ... , 



and the subset is complete in L 2 ^ x ' 



jp S to ^umplete in L 2 *yux, j. 

As usual, if {ipp} is the adjoint basis of eigenfunctions of the adjoint operator 

(4.6) B = -A 2 + iyV + f I in L 2 p (R N ), with p=i 

with the same spectrum (14.41) . the bi-orthonormality condition holds in L 2 (M N ): 

(4.7) (ijjft, ip*) = 5p V for any fj,, v. 

4.2. Inner expansion. Thus, in the Inner Region characterized by compact subsets in 
the similarity variable y, we assume a centre or a stable subspace behaviour as r — > +oo 
for the linearized operator B* + /: 

centre: Y (y , r) = a(r)ty(y) + w 1 (\p = -1, \(3\ = 4), 

' ' stable: Y(y, r) = -Ce x ^(y) + w 1 (\p < -1, \(3\ > 4). 

For the centre subspace behaviour in (14.81) . substituting the eigenfunctions expansion into 
equation (14.21) yields the following coefficient: 

(4.9) a = a 2 7o + where 7o = c ((^) 2 , ijip) => a(r) = + ... . 
Note that for the matching purposes, we have to assume that (see details in 

(4.10) if i/j*(0) > 0, then 7o > and C > 0. 
Actually, for k = 4 and iV = 1, it is calculated explicitly that 

(4.11) 7o = -c 136V6 < 0, 
so that the centre manifold patterns with the positive eigenfunction 

(4.12) n{y) = -h ^ 4 + 24 ) 



■>(, 



correspond to solutions that blow-up on finite interfaces; see [iQl § 3]. A full justification 
of such a behaviour can be done along the lines of classic invariant manifold theory (see 
e.g. [60]), though can be very difficult. 

Actually, we can construct more general asymptotics by taking an arbitrary linear 
combination of eigenfunctions from the centre subspace. Overall, the whole variety of 
such asymptotics is characterized as follows: 

f -±- + ... for \3\ = 4, 

(4.13) Y(y,T) = a(T)\yf X (v) + ..., where a(r) = - + for ™ > ^ 

In general, here, x{f) > is an arbitrary smooth function on the sphere S^ -1 , where its 
positivity is induced by matching issues to be revealed below. 

4.3. Outer region: matching. We follow [12], where it is shown that the asymptotics 
(14. 8 p admit matching with the Outer Region, being a Hamilton-Jacobi (H-J) one. More 
precisely, in the centre case with \/3\ =4, according to (I4.8p . (14. 9p . we introduce the outer 
variable and obtain from (13. 2p the following perturbed H-J equation: 

(4.14) £=7^4 =► u^-i^Vv-^u + H^v + iQe-Vv-A 2 ^. 

Passing to the limit as r — > +oo in such singularly perturbed PDEs is not easy at all 
even in the second-order case (see a number of various applications in [39]). Though 
currently not rigorously (this looks being completely illusive), we assume stabilization to 
the stationary solutions /(£) satisfying the unperturbed H-J equation: 

(4.15) V/-^/+|/rV = in R N . 

This is solved via characteristics, where we have to choose the solution satisfying (I4.13P : 

(4.16) /(0 = /*-^|?| 4 xM + - as =► /(0 = /*(l + c*|e| 4 xM)"^ T , 
where c* = (p — l) p /( p_1 ). Since 7 < according to (14.111) . the resulting profile 

satisfies /(£) > /* and blows up on the surface {c* l^^xC^) = ~~ !}• Note that this actual 
nonexistence of a bounded centre subspace pattern plus the known unstable eigenspace 
of B* + I in the linearized equation (14.21) somehow reflect existence of two self-similar 
solutions fo(y) and fi(y) as "nonlinear eigenfunctions"; see [40J. 

Thus, these centre subspace patterns are not bounded and should be excluded from the 
consideration. On the other hand, for the 2mth-order PDE (11.171) with odd m = 3, 5, ... , 
we have 70 > and then (I4.16P can represent standard blow-up patterns, [i~9] . 

Similarly, for the stable behaviour for > 4 in (14.81) . we use the following change: 

4-l/3| |/3|-4 

(4.17) t_ = e m T y v T = - pi ^-Vv - ^v + lv^v-e T A 2 v. 

Passage to the limit r — > +00 and matching with Inner Region are analogous and lead to 
truly existent blow-up patterns for any m > 2 in (11.171) . [19]. 

Overall, according to matching conditions (I4.13P and (14.171) . the whole set of possi- 
ble blow-up patterns of Type II(Her) is composed from a countable set for \j3\ = 4 (m 
odd), 5, 6, ... of continuous (uncountable) families induced by smooth functions x on S 7 ^ 1 - 
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5. Type Il(sing): LINEARIZATION ABOUT THE SSS AND MATCHING 

The idea of such Type II blow-up patterns for the RDE-2 (11.71) is due to Herrero- 
Velazquez [32], where a justification of existence was achieved (see [HS] for extra details). 
We apply this method to the RDE-4 (11 .11) and, by the same reasons, we are not obliged 
to concentrate on a proof. Thus, instead of the linearization (14.21) about the constant 
equilibrium, we perform it about a singular one. 

5.1. Singular stationary solution (SSS). Consider the stationary equation 
(5.1) - A 2 U +\U\ p - 1 U = in R^UO}. 

The explicit radial SSS has the standard scaling invariant form 

(52) U(y) = C,\y\-», where ft = ^, C* = , 

and D = fi(fi + 2)[(ti + l)(fj, + 3) + {N-l)(N-5-2fi)]. 
It follows that such SSS exists, i.e., D > 0, in the following parameter ranges: 

(5.3) p >£_, iv>4, or p<g JV>2. 

5.2. Linearization in Inner Region: discrete spectrum by Hardy— Rellich in- 
equality. Thus, we perform linearization in (I3.2p about the SSS: 

(5.4) v = U + Y v T = B*y + D(y), 
where, as usual, T)(Y) is a quadratic perturbation as Y — > and 

(5.5) B* = H*- \y V-^J and H* = -A 2 + $ I. 

Similar to Lemma 14.11 the operator B* at infinity admits a proper functional setting in 
the same metric of L 2 * . However, it is also singular at the origin y = 0, where its setting 
depends on the principal part H*. 

Proposition 5.1. The symmetric operator H* admits a Friedrich's self-adjoint extension 
with the domain Hq(B\), discrete spectrum, and compact resolvent in L 2 (Bi), where Bi C 
M, N is the unit ball, iff 

(5.6) pD < ch = 

Proof. Indeed, (I5.6P is just a corollary of the classic Hardy-Rellich-tjpe inequality 

(5.7) m ^Jj^<J\M 2 for ueH 2 ( Bl ), 

where the constant is sharp. For compact embedding of the corresponding spaces, see 
Maz'ya [621 p. 65, etc.]. □ 



4 This was derived by Rellich already in 1954; see [13] and [55] for further references and full history. 
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p=2: the N-range of applicability of Hardy's inequality 
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p=3: the range of applicability of Hardy's inequality 




(a) p = 2 (b) p = 3 

FIGURE 25. Checking inequality (pTHl) : p = 2 (a) and p = 3 (b). 



(5.8) 



The necessary inequality (15. 6p takes the form 

x[(i + Fi)(3 + ^r) + (iV-i)(iV-5-^)] >0 

and does not admit an easy analytic analysis. In Figure I2"S*| numerics show that 

(5.9) ([E2D holds for TV > 24 if p = 2, and iV > 19 if p = 3. 
In particular, checking (15.81) at p = +oo yields the inequality: 

(5.10) G^N) = lN{N w i)]2 - 8[3 + (N — 1)(N — 5)] > 0. 
If this is true, then ( 15.81) holds for all p> 1, so: 

Proposition 5.2. For any N > 13, there exists a pn(N) > 1 such that 

(5.11) (123) holds for all p > Ph(N), 

and hence the operator B* in (I5.5|) has a discrete spectrum in L 2 , i 



5.3. Inner Region I. Thus, we assume that, under certain conditions, (15.111) holds and 
cr(B*) = {Xk} is discrete, with the eigenfunctions {ipp, = k}. Furthermore, it is also 
convenient to assume that the spectrum is (at least partially) real. To justify such an 
assumption for this non-self-adjoint operator, we rewrite (15.51) in the form 



(5.12) 



B* 



B* + |4r / /, where c = pD 

\yr p-1 ' L 



and B* is the previous operator ( 14. 3 j) with the real spectrum shown in Lemma 14.11 (ac- 
tually, this means that B* admits a natural self-adjoint representation in the space l 2 p « of 
sequences, where it is also sectorial, [25]). Therefore, the real spectrum of (15 . 12[) can be 
obtained by branching-perturbation theory (see Kato [51]) from that {\p = — | — k = 
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\(3\ > 0} of B* 1 at c = 0. Next, the branch must be extended to c = pD, which is 

also a difficult mathematical problem; see [2H § 6] for some extra details, which are not 
necessary here in such a formal blow-up analysis. 

Thus, we fix a certain exponentially decaying pattern in Inner Region I: 

(5.13) Y(y, r) = Ce X P T fy(y) + ... as r -> +oo (Ag < 0). 

If there exists A = G cr(B*), the expansion will mimic that in (14.81) for the centre 
subspace case. Note that (I5.13P includes all the non-radial linearized blow-up patterns. 

5.4. Matching with Inner Region II close to the origin. In order to match (15.131) 
with a smooth bounded flow close to y — 0, which we call Inner Region II, one needs the 
behaviour of the eigenfunction i>%(y) as y — > 0. To get this, without loss of generality, we 
assume the radial geometry. Then, the principal operator in the eigenvalue problem 

(5.14) H*V>* + ••• = A^* as y -> 
yields the following characteristic polynomial (see |23j): 

(5.15) r(y) = \yy + ... H c (n) = - 7 (7-2)(7 + iV-2)(7 + A r -4) + c = 0. 
Consider the most interesting critical and extremal case 

(5.16) c = pD = c Bi =^^- =► Hc ( 7 ) = -b + ^[f + (N-A) 7 -^]. 
There exists the double root 7^2 = — < 0, which generates two L 2 -behaviours: 

(5.17) Vi(y) = \y\-—\n\y\(l + o(l)) and ^(y) = 1 2/ 1 ^ (1 + o(l)) as 

Note that ifp- approximations of establish that ch is the best constant in (15.6jl . Other 
two roots of the characteristic equation in (15.161) are 

(5.18) 73,4 = \ [4 - N T v/(iV-4) 2 + iV 2 ] , 
where 73 < 71,2 < and 74 > corresponds to L 2 -solutions. We have 

(5.19) r 3 (y) = \y\^(l + o(l))?L 2 , 

so that in L 2 the deficiency indices of B are (3,3) and cannot be equal to (4,4). Unlike 
the second-order case, the straightforward conclusion on the discreteness of the spectrum 
in the case (4, 4) [TQl p. 90] does not apply, so Friedrich's extension of H* is constructed 
by other arguments [23J and include settings, where two most singular behaviour in (15.171) 
and (15.191) are excluded. 

Overall, this gives the following behaviour of the proper eigenf unctions at the origin: 

(5.20) ^%{v) — — ^M - "^ + ••• as V — * (vp > are normalization constants). 

This allows to detect the rate of blow-up of such patterns by estimating the maximal 
value of the expansion near the origin: 

(5.21) vp(y, t) = C*\y\~~ — vpCe^ l3T \y\~ E z ± + ... as y — > and r — ► +00, 
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where we observe the natural condition of matching: 

(5.22) vpC > 0. 

Calculating the absolute maximum in y of the function on the right-hand side of (15.221) 
(this is a standard and justified trick in some R D problems; see e.g., [8]) yields an 
exponential divergence: 

(5.23) |M-,r)|| 0O = ^e^ + ..., where pp = {N _f { ; l _ ps) > {p>Ps), 

where dp > are some constants. Note that, depending on the spectrum {A^ < 0}, (15.231) 
can determines a countable set of various Type II blow-up asymptotics. 

Let us define more clearly the necessary matching procedure. In a standard manner, 
we return to the original rescaled equation (13.21) and perform the rescaling in Region II 
according to (15.231) : 

(5.24) v(y,r) = e^w(C,s), £ = ^y, »p = a = e^" 1 ^. 
Then w solves the following exponentially perturbed uniformly parabolic equation: 

(5.25) w s = -A 2 w + \w\^w - i [(J + n p )Z .V w+ (-L l + Pp ) w ] . 

As above, we arrive at a stabilization problem to a bounded stationary solution, which is 
widely used in blow-up applications (see examples in [39J). In general, once the uniform 
boundedness of the orbit {w(s), s > 0} is established, the passage to the limit in (15.251) 
as s — > +oo is a standard issue of asymptotic parabolic theory. 

Our blow-up patterns correspond to the stabilization uniformly on compact subsets: 

(5.26) w(£,s)^W(£), s^+oo, where - A 2 W +\W\ P ~ 1 W = 0, £ G M. N , W(0) = dp, 
for all admissible \/3\ — 0, 1, 2, ... . We next discuss a crucial issue on such a matching. 

5.5. Matching: on necessary structure of global bounded stationary solutions. 

There are two issues associated with the stationary problem ( 15.26)) . 

1. Firstly and elementary, one can see that, bearing in mind the matching of Regions I 
and II, the bounded stationary solutions W(£) defined by (15.261) must be positive and non- 
oscillatory as £ — » oo. Otherwise, such a matching with positive SSS U(£) is impossible. 
There exists a definite negative result in the subcritical Sobolev range (there is a diverse 
literature on this popular nowadays subject, so we refer to a recent paper [15] as a guide): 

(5.27) a solution W > of (I5.26P is nonexistent for p e (l,ps). 

Actually, this means that all the entire (i.e., without singularities) solutions of (I5.26P are 
oscillatory, as Figure |2"61 shows for p = 2, iV = 3. Note that we are restricted by (15.61) . 

2. Secondly and fortunately, existence of such positive solutions W(£) is well established 
already [H p. 908]: 

(5.28) for p > p s , for any dp > 0, there exists a unique positive solution W(£). 

Here we exclude the critical case p = ps, where exact positive solutions exist to be used 
in Section [61 As a numerical illustration, Figure [27] shows two such results for iV = 13 
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FIGURE 26. Shooting an oscillatory solution of (I5.26P for p = 2 and N = 3. 

(a), where the dotted line denotes the explicit solution for N = 12. It is clearly seen that 
W(£) for N = 13 lies above this, so remains positive. In (b), we show the positivity of 
the solution W for N = 24, where by (|5.9p . the spectrum is guaranteed to be discrete. 

3. Of course, the above is not sufficient for matching of Inner Regions I and II to get 
a blow-up pattern. More importantly, we have the following: 

Proposition 5.3. The entire solutions W(£) of the radial ODE (I5.26P are not oscillatory 
as y — > +oo about the SSS (15. iff (15. iil) holds, and then: 

(5.29) p > '■ W(0 has at most finite intersections with U(£) on £ G (0, +oo). 

Proof. It suffices to observe that, as customary, the oscillatory behaviour as y —>■ +oo is 
governed by the linearized operator therein, which is (15.51) (the limit (15.311) below justifies 
the linearization). Hence, in the critical Hardy case, the characteristic polynomial (15.161) 
has real roots only (actually, all of them, and this is quite a general property [23, 24J), 
and obviously the same holds in the sub critical range pD < ch, meaning that W(£) is not 
oscillatory about U(£) as £ — ► +oo. Clearly, if pD > ch, (15.151) and (15.161) imply existence 
of a proper root 7 G C with a not that large negative real part. □ 

Thus, we have concluded that, for the present problem: 

(5.30) discrete spectrum and non-oscillation occur in the same range p > pn(N). 

Indeed, this has some natural roots in general spectral theory of ordinary differential 
operators. For instance, for second-order singular operators, the non-oscillating behaviour 
at singular endpoints always imply existence of a self-adjoint extension in L 2 with a 
discrete spectrum; see Lemma 3.1.1 in [58, p. 74]. For higher-order symmetric operators 
[70] . such a universal conclusion is not that clear, though is easily observed in particular 
problems related to simpler homogeneous operators for Hardy's inequalities as in [23| 124"]. 
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Proposition 15.31 for p > pn(N) was proved in [42, p. 909], where other important 
properties of entire solutions W(£) of (I5.26P have been established. So we do not need to 
mention them here in detail and will use the following only (see also [7j): 

(5.31) (£±i)F ± r>m^i as f_> +00 . 

However, a number of problems concerning (15.2ft remain open. For instance, proving that 
(cf. Open Problem 3 in [42, p. 915] on ordering of the family {W(£), dp > 0}) 

(5.32) for p > pn(N), W(£) does not intersect U(Q. 

Note that in view of inevitable using shooting techniques, the property (I5.32p is very 
difficult to check numerically. 

Fortunately, as a standard topology suggests, solving the open problem (I5.32|) is not nec- 
essary for the validity of the matching of Inner Regions I and II, since the non-oscillating 
of W(£) as £ — > +oo is in principal demand (one can see that existence of a finite number 
of intersections cannot spoil matching). Thus, we conclude that: 

(5.33) for p > ph(N), matching of two flows (15.131) and (I5.26P is plausible, 

though a huge mathematical work is necessary to prove this (the author still believes that 
this can be done in a reasonably finite period of time, but its scale can be beyond any 
expectation). 

5.6. On new blow-up similarity solutions in the oscillatory range p < pn- Thus, 
( 15.161) clearly shows that for p < pn the solutions W(£) are oscillatory about the SSS £/(£). 
Such topology (as in the second-order case, see [38] and later publications) suggests that 
in this subcritical Hardy range there may be a sequence of similarity profiles satisfying 
(12.171) and exhibiting arbitrary finite oscillations about U(£) for sufficiently small radial 
£ > 0. Such self-similar blow-up profiles concentrated in a neighbourhood of the unstable 
singular equilibrium U (above U, a.a. solutions must blow-up), are expected to be also 
highly unstable, at least in comparison with the previous profiles /o and f\ studied in 
Section [2J Therefore, we ignore such new families (possibly countable depending on 
parameter ranges) of the s-s blow-up. 

5.7. On related non-radial blow-up patterns. These can be predicted in a various 
ways. Firstly, one can start with a non-radial SSS solving the elliptic equation (15. ip . but 
surely such ones are unknown. Secondly, under the condition (15.61) . a non-radial eigen- 
function ip%{y) (e.g., corresponding to an "angular" swirl obtained by angular separation 

of variables) of B* can be taken into account. Then matching will assume using non-radial 
entire solutions of (I5.13p . which then deserves further study. 
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Radial stationary problem -A^V + |Vf~'V=0: p=2, N=13,C,=840 



Radial stationary problem -A 2 V + |V| P 'v=0: p=2, N=24 (Hardy range) 
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(a) N = 13 > 12 = N 2 




(b) N = 24 



FIGURE 27. Checking the positivity of the solutions of (15.261) for p = 2: iV = 13 
(a) and iV = 24 (b). 



6. Type I(LN): NON SELF-SIMILAR BLOW-UP EVOLUTION ON A MANIFOLD OF 
GENERALIZED LOEWNER— NlRENBERG STATIONARY SOLUTIONS 

6.1. Classic Loewner— Nirenberg (L— N) conformally invariant exact solutions. 

These are classic solutions obtained in Loewner-Nirenberg [59] in 1974 for the second- 
order elliptic equation 



(6.1) 



AW + W p = in 



W(0) = d>0, for p = Ps 



N+2 
N-2 



{N>2), 



which are invariant under conformal and projective transformations (symmetries of (16.11) 
were earlier studied by Ibragimov in 1968 [51]). These solutions are given by 



(6.2) 



Wb(0 = d 



N(N-2) 



N-2 
2 



> in 



N(N~2)+d i /( N ' 2 )\^\ 2 

and exhibit a number of uniqueness and other exceptional properties of the equation (16. ip . 

6.2. Generalized L— N solutions for the biharmonic equation. For the critical bi- 
harmonic counterpart of (16.11) 



(6.3) 



-A 2 W +\W\ p - 1 W = 



in 



W{0) =d>0; p = p s 



N+4 
N-4 



(N>4), 



the corresponding exact solutions are known from the 1980s at least, which we call the 

generalized L-N ones: 



(6.4) 



W (0 = d 



v /(iV+2)(jV 2 -4)(jV-4) 



> in 



. v /(Af+2)(AT 2 -4)(AT-4)+ ( i 4 /(iv-4) 

The earliest references to the exact expressions (16. 4p we have found are [33l p. 1057] in 
1985 and [72l [82] in 1992, where in the latter one important properties of Wq have been 
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proved (see also [45 J for further references). Note that, for the 2mth-order polyharmonic 
extension, the corresponding positive entire solutions look similarly: 



-(-A) m W + {W^W = 0, p = f±g =}► W (0 = d 

where N > 2m and B m = [jv^lm+i))!! • ^ ee Svirshchevskii in 1993 [81] (in a preprint, the 
solutions were published as earlier as in 1989 [80j), and more related exact solutions of 
other critical elliptic PDEs (e.g., with a p-Laplacian) and extra references in [321 § 5]. 

6.3. Formal construction of Type II(LN) blow-up patterns for p = p s . Let v(y, r) 
be the rescaled solution of (I3.2p in, say, radial geometry at the moment. Let us assume 
that v(y,r) behaves for t ^ 1 being close to the stationary manifold composed of the 
explicit equilibria (16.41) . i.e., for some unknown function ip(r) — > +oo as r — > +oo: 

(6.5) v(y,T)=(p(T)W {tp x Tr(r)y)+... 

on the corresponding shrinking compact subsets in the new variable ( = Lp^(r)y. It 
then follows that, on the solutions (16.41) in terms of the original rescaled variable y (cf. 
computations in [UJ p. 2963]; our notations have been slightly changed) 

(6.6) \v(y,T)r 1 v(y,r)^^d(y) as r - +oo 

in the sense of distributions, where > is some constant. Therefore, on this manifold 
of solutions, the rescaled equation (13.21) takes asymptotically the form 

(6.7) v T = A(v) = -A 2 v-ly-Vv-^v + ^5(y) + ... for r » 1. 
According to Lemma 14.11 we are looking for Type II patterns of the form 

(6.8) Vf,{y, t) = c p {r)i)*p{y) + ... c p = -a p c p + + ... , 

where ap = 2 I /3 I+ Af ~ 4 > o and hp = eNipp(0)- Simple particular "resonance" solutions 
correspond to an exponential divergence: 

(6.9) (p(r) =e Q " r + ... and cp(r) = hpre' a ^ T + ... for r > 1, > 0. 

Bearing in mind the scaling in (16.51) . this yields a countable family of distinct compli- 
cated blow-up structures, where most of them are not radially symmetric. To reveal the 
actual space-time and changing sign structures of such Type II patterns, special match- 
ing procedures apply. In [TJ], this analysis has been performed in the radial geometry 
for (11.71) . though still no rigorous justification of the existence of such blow-up scenarios 
is available. Thus, the first Fourier coefficient in (16.81) implies a complicated structure 
of the pattern around the formed Dirac's 8(y) according to (16.61) . However, since these 
expansions are given by generalized Hermite polynomials {ipp} this matching is expected 
not to impose more difficulties as those similar in Section HI In any case, more matching 
details for the much harder PDE (II. ip seem then excessive here. 
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Appendix A: On universality of the open L p =>• L°° problem in PDE theory 

The Millennium Prize Problem, posed specially for the NSEs (|1.18p . is, in a loose sense, "non- 
unique" , since similar open regularity problems (or not that lighter significance) occur for many 
evolution PDEs of various types. We list a few of them, where the difficult open mathematical 
aspects of global existence and/or blow-up are associated with the following factors: 

(i) supercritical Sobolev parameter range of the principal operator (hence, standard or very 
enhanced embedding-interpolation techniques fails), and, in fact, as a corollary, 

(ii) multi-dimensional space x £ R , with A^ > 3, at least (this leaves a lot of room for 
constructing various L°° blow-up patterns via self-similarity, angular swirl, axis precessions, 
linearization, matching, etc.). 

We now list those PDEs, where we give a few recent basic references to feel the subject. 

(I) Supercritical defocusing nonlinear Schrodinger equation^ (NLSE) (see [B3]-[II5], [851 [75]) 

(A.l) -iu t = Au- |u|P-y with p > p s (2) = $±§ ( N > 3 ); 

(II) 2mth-order supercritical semilinear heat equation with absorption (m = 1 is covered by the 
MP; see [41] and [6], where the result in § 4 for p > ps(2m) applies to small solutions only): 

(A.2) ut = -{-A) m u - \u\ p ~ l u, with p > p s (2m) = g±|a (N > 2m, m> 2); 

(III) The semilinear supercritical wave equations (see [52|, [86] , as most recent guides) 

(A.3) u tt = Au - MP-y with p > p s (2) = (A > 3). 

Possibly, here the Maximum Principle kind arguments associated with the single Laplacian A 
can still play a role; then it is to be replaced by —A 2 ; see below. 

One can add to those "supercritical" PDEs some others of a different structure such as the 
Kuramoto-Sivashinsky equations for I = 1,2,.. 



(A.4) u t = -(-A) 2l u + (-A) l u + ± E( k )dkD Xk (\un |d| = 1, p> Po = l + ^ 21 . 

Here, po is not the Sobolev critical exponent, though precisely for p > po, I? ^ L°° by blow-up 
scaling, |34[ § 5]. On the other hand, a more exotic applied models exhibit similar fundamental 
difficulties such as the following nonlinear dispersion equation (see |26[ [35] for references and 
some details) 

(A.5) u t = -D Xl [(-A) m u}-D Xl (\u\P~ 1 u), with p > p s (2m) = #±|2 

In view of the conservation properties for the models (|A.4p and (|A.5p . these, though being local, 
can be more adequate to the nonlocal NSEs (|1.18p . than the others above. 

In most of the cases, the operator on the right-hand sides satisfying for u G C^°(R N ) 

(A.6) A{u) = -(-A) m u-\u\P- l u => {A(u),u) = - f\D m u\ 2 - J\u\p +1 <0, 

is indeed coercive and monotone in the metric of L 2 (M. N ), which always helps for global existence- 
uniqueness of sufficiently smooth solutions of these evolution PDEs. For the NLS (jA.ip . this gives 
a stronger conservation laws than for the focusing equation with the "source-like" term = |u| p_1 n. 



'The author would like to thank I.V. Kamotski, who first attracted his attention to this problem. 

40 



Evidently, replacing A in (|A.ip and (|A.3p by — (— A) m , m>2 moves the supercritical range to 
that in (IA.2j) . On the other hand, introducing quasilinear differential operators — {—A) m \u\ a u 
with a > moves the critical exponent to ps(2m, a) = (cr + 1) jv-2m ' Similar supercritical PDEs 
can contain 2mth-order p-Laplacian operators, such as the one for m = 2, with a > 0, 

(A.7) A(u) = -A(\Au\ a Au)-\u\P~ 1 u, {A(u),u) = - f \Au\ a+2 - f \u\p +1 < 0. 

However, the lack of embedding-interpolation techniques to get L°°-bounds, which can be 
expressed as the lack of compact Sobolev embedding of the corresponding spaces for bounded 
domains Q C (this analogy is not straightforward and is used as a certain illustration only) 

(A.8) H m (Q) <£_ LP +1 (Q) for p > p s {2m), 

actually presents the core of the problem: it is not clear how and when bounded solutions can 
attain in a finite blow-up time a "singular blow-up component" in L°°. For the operator in 
(IA.7j) . a similar supercritical demand reads 

(A.9) W*+ 3 (n) t H*i(n) for p> Ps (A,a)= ( ^%f 2) , ^> 2(<r + 2). 

In the given supercritical Sobolev ranges, finite mass/energy blow-up patterns for (|A.ip - (|A.5p 
are unknown, as well as global existence of arbitrary (non-small) solutions. 
It is curious that for the NSEs with the same absorption mechanism as above, 

(A.10) u t + (u- V)u = -Vp + Au- |u| p - 1 u, divu = in R 3 xR + , 

by the same reasons and similar to (|A.2p , the global existence of smooth solutions is guaranteed 
[3] in the subcritical Sobolev range only: for 

(A. 11) p < 5 = |£z§|jv_3 (and p > | by another natural reason). 

We thus claim that, even for the PDEs with local nonlinearities (|A.ip - (|A.3P (and similar 
higher-order others), the study of the admissible types of possible blow-up patterns can represent 
an important and constructive problem, with the results that can be key also for the non-local 
parabolic flows such as (|1.18p . (|A.10p . etc. Moreover, it seems reasonable first to clarify the blow- 
up origins in some of looking similar and simpler (hopefully, yes, since (I1.18P is both nonlocal 
and vector-valued unlike the others) local supercritical PDEs, and next to extend the approaches 
to the non-local NSEs (j!.18f) : though, obviously, the former ones are not that attractive and, 
unfortunately, are not related to "millennium" issues (however, many PDE experts very well 
recognize how important these are for general PDE theory). 
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